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Aerodynamic Maneuvering Hypersonic Flight Mechanics 


1. Introduction 

The emergence of current high-interest missions involving aeromaneuvering 
hypersonic flight has given rise to the corresponding need for preliminary 
design and performance analyses of such vehicles. This need in turn has 
motivated efforts to develop simplified analytical and computational methods for 
parametric analysis of maneuvering hypersonic fiight under conditions 
appropriate to the mission involved. The missions of interest include those of 
aircraft, e.g., the Orient Express, and hybrid or •'spaceplane" vehicles, e.g., 
the singie-stage-to-orbit or transatmospheric Vehicle (TAV), and the 
aeroassisted spacecraft such as the Aeroassisted Orbital Transfer Vehicle 
(AOTV). These missions feature aeromaneuvering flight at high altitudes and 
high Mach number, and are characterized by protracted global-scale maneuvers. 
Under these conditions, the trajectory may be significantly affected by 
rotating-Earth and Earth-curvature effects in addition to the primary balance of 
gravity, thrust, and aerodynamic forces. In addition, the characteristic 
trajectory parameters as well as the associated coordinate system in which the 
trajectory is described require evaluation since aircraft related parameters and 
frames of reference are typcally different from those of oroiting spacecraft and 
reentry trajectories. 

The purpose of this study was to develop simplified analytical methods for 
parametric analysis of hypersonic maneuvering flight. The effort included a 
review of different formulations of the general equations of motion, their 
associated coordinate frames, various simplifications of the equations, and 
previously achieved analytical solutions. This study sought to both extend 
previous solution methods and to develop new ones. In addition, evaluation of 
the literature and developing a systematic perspective on the knowledge it 
represents proved to be a major portion of the effort. 

2. Background 


2.1 Literature Categorization 

A review of the hypersonic maneuvering-f light literature involves 
categorization into planar/nonplanar studies, studies including/excluding 
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propulsive thrust, and studies using an Earth-fixed wind axes coordinate frame 
or Earth-centered spherical coordinate frame or orbital axes/elements frame. 
Also, the varied assumptions of maneuver constants and equation transformations 
employed in the previous studies are of interest. Excluded from this review are 
earlier studies involving ballistic reentry vehicles and studies focussing on 
the guidande/navigation/control optimization of hypersonic maneuvering flight. 
References 1-3 are examples of general studies presenting full derivation of 
the equations of motion with both analytical and computational solutions. 
Shkadov, et all treat prqpelled aircraft as well as thrustless maneuvering 
descent from orbit, both with Coriolis effects excluded. Loh2 treats the 
various cases of thrustless lifting and nonlifting reentry, again with Coriolis 
effects excluded from analytical solutions. Vinh, et al3 derive the full 
equations of motion, and then exclude from analytic solutions Coriolis effects 
as well as thrusting during reentry. The scarce literature dealing with 
thrusting maneuvering hypersonic flight is found in References 4 through 13. 
References 4-10 focus on the nonplanar (turning) "aerocruise" maneuver where 
thrust continually cancels drag. References 11 - 13 deal with simplified 
analytical solutions for planar maneuvering Transatmospheric Vehicle (TAV) 
flight. The prolific nonthrusting or "aerogliding" literature includes 
References 13 through 27 and 36 for nonplanar flight and References 20 through 
34 for planar flight. Reference 35 provides a general theoretical formulation 
applicable to all kinds of missions. 

By far the majority of the literature employs a wind axis coordinate frame 
in which the equations of motion are expressed. As exceptions, equations are 
expressed in a spherical coordinate frame in References 5, 8, 15 and 29, and in 
an orbital element frame in References 9 and 10; the frames employed are not 
obvious in References 12 and 24. Finally, pure computational solution 
techniques are used in References 6, 7, 12, 19 through 24, 29, 30, 33 through 
36. All other references to some extent provide analytic solutions, some with 
computational solutions for comparison. None of the references employing wind 
axis coordinates and providing analytical solutions include Coriolis effects on 
the trajectory and vehicle performance. Further, analytic-solution references 
show a bewildering variety of parameters chosen as the independent variable of 
differentiation in the transformed differential equations of motion that 
facilitate analytical solutions, variables are used such as flight path 
angle y, time t, and distance along the trajectory s. These analytic references 
also show a large variety of assumptions on constant parameters during 
maneuvering flight: References 4 and 7 hold velocity V and altitude (and thus 


density) constant during the turn, References 4, 6 and 7 bank angles', and 
References 4, 7, 11 and 13 dynamic pressure q. 


3. General Equations of Motion - Different Formulations 
3.1 Earth-f ixed Geographical Frame 


Vinh, et all provide a general formulation for the nonplanar equations of 
motion in an Earth-fixed rotating frame where thrust and aerodynamic forces 
(lift and drag) are included. The reference plane is the equatorial plane, so 
that vehicle position is described in latitude/longitude/altitude coordinates 
(geographical coordinates). These equations comprise a determinate set of three 
dynamical and three kinematical equations in six unknowns as follows: 


at _ a -v Cl) 4 * Cos (s 'Vfr c OS <1> — 

At ** 0 J 

V = Tj/wf.+L c tte'—fa— — t-t 2<*>Vcot 6co* <t> + aj 9 t- 

C* U * 

V 4z “ 5125 - JL’caiTfco-ty +• &l»<PcAS<p cos V 


a V c tn Y 
dt 

d & VcwYtw V' 


t; d, c * «/ ■ft" 

£ = 'ftlrus'h Qnj/e lij A+ 

kvi = »*CLSS 

C~= hall C^aiHk.) atyle 
(*> - h> ' t«.+/6n hafe. 

CV - s^e«4 

O«k*ou/*.S «^<c C t is up) 

V— OUkwi/R Uetldinj 0^v***\ EVsi~) 

W a 4«si'ho n 

0-w la+i'VuJe 
((is lovijI^uJe 


The equations assume a spherical, rotating Earth, and out-of-plane force is 
generated by banking the lift force. In the discussion of these equations, the 
various contributing "acceleration" terms or effects are referred to as lift 
(L/m), drag (D/m), thrust (T/m), gravity (g), centrifugal (v2/r), Coriolis 
<2toV), and Earth-centrifugal (ku 2) . The last two are "Earth-rotation" 
effects, while the centrifugal term v2/r is due to flight path curvature. 
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Table 1 estimates the order of magnitude (in g units) of each of these 
contributing terms for a variety of missions. It becomes obvious that although 
the Earth centrifugal terms (kj£) may be neglected for all missions, the 
Coriolis term should be retained for high-altitude, hypersonic global-scale 
missions such as the TAV and the AOTV. This term is not retained in any of the 
earlier studies reporting analytic solutions to the equations. The development 
of corresponding analytical solutions where the Coriolis term is retained has 
been a chief focus of the present work. It is of particular interest to compare 
accelerations in g units as a function of tlach number for (gravity + 
centrifugal), Coriolis, and Earth-centrifugal effects (Table 2) It is seen that 
the Earth-centrifugal effect is always neglectable, but Coriolis is 10% or more 
of the sum (gravity + centrifugal) at tlach numbers of 15 and above (about 5% at 
li — 8) . 

Finally, Equations 1-3 degenerate to those employed in much of the 
planar-flight analytic-solution literature (References 11 - 13) when all the 
Earth-rotation effects are neglected (and E, O' = 0): 

( 4 ) 4^. — — - — -*\s /*"£ 

^ At *vs ^ d 

dt ♦V' J 3*" 

The Vinh formulation involves Earth-fixed latitude/longitude, i.e., 
geographical coordinates, and is a natural formulation for Earth-relative motion 
in nonmaneuvering reentry. The equatorial plane is the reference plane, making 
downrange/crossrange maneuvers very difficult to express and visualize by means 
of this formulation. 

3.2 Orbital Frame 

CervisilO presents a succinct summary of the equations for motion assuming 
a spherical nonrotating Earth and expressed in terms of the variation of orbital 
elements caused by aerodynamic forces. These equations are coupled, 
first-order, nonlinear differential equations in terms of time derivatives of 
the orbital elements: 
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In the absence of lift, drag, and thrust forces, these equations describe a 
conic orbit. Analytic and computational solutions built on this framework have 
largely been restricted to aerocruise and gliding plane-change (inclination 
change) maneuvers whose application of interest is to AOTV-like missions. 

Clearly, this formulation expresses the modification of orbital motion 
produced by transient aerodynamic forces. It is a natural formulation for 
spacecraft maneuvers involving skip- trajectory passage through the atmosphere. 

It uses as the reference plane the original orbital plane, which is 
Earth-centered (and equivalent to a great-circle ground track). Thus, 
maneuver-produced downrange/crossrange deviations are more easily expressed with 
this formulation than with the previous one. However, the motion is inertially 
referenced, so that Earth-relative motion is not easily expressed or visualized 
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3.3 Crossrange/Downrange Earth-Fixed Formulation 

Ikawa35 has provided a useful formulation of the equations of motion that 
expresses Earth-relative motion in coordinates that use as the reference plane 
the initial orbital plane (or heading plane for atmospheric vehicles). This 
formulation leads naturally to downrange/crossrange motion solutions and 
visualization. As such, it is probably the most useful formulation for 
aeromaneuvering missions. The formulation describes orbital (conic) motion in 
the absence of thrust and aerodynamic forces. It was implemented on an IBI1 PC 
host microcomputer, and Ikawa describes accurate solutions for several 
distinctly different missions: 1) orbital motion and ground track for two 
different oroits including a 24-hour, 600 inclination orbit; 2) reentry 
missions including a synergetic orbit plane-change maneuver, and a more general 
mission of aeorbit-to-lano from elliptic LEO using aeromaneuvering; and 3) a 
cross-country subsonic commercial aircraft flight. 

The general equations of motion developed in this formulation are as 
follows : 
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4. Solutions to Forms of the General Equations 

4.1 General Solutions 

Trajectory solutions have been achieved computationally in the past for all 
types of missions - orbital motion and maneuvers, aeromaneuvering skip 
trajectories (orbit-to-orbit) , nonmaneuvering and maneuvering or bit- to-ground 
(reentry), single and multi-stage ascent to orbit (maneuvering and 
nonmaneuvering), and hypersonic maneuvering aircraft flight. Hypersonic 
aircraft, ascent, and orbit-to-ground (reentry) simulations typically utilize 
the wind-axis, Earth-relative geographical coordinates formulation, e.g., 
Equations 1-3. Orbit-to-orbit simulations either use this same formulation or 
the perturbation-in-orbital-elements formulation, e.g., Equations 6 - 11. 

It would appear that all maneuverable misssions, except perhaps 
orbit-to-orbit, would more easily be visualized and studied parametrically when 
expressed in downrange/crossrange motion coordinates, and this is the utility 
afforded by the formulation of Ikawa.35 

Analytic solutions to simplified forms of the general equations have also 
been developed in many studies over the years. The most prevalent (and most 
relevant) solutions are those developed for forms of the planar wind-axes 
equations under various flight conditions. To describe and summarize these 
literature results would require a very lengthy discussion. 

Instead, the few results achieved in the present study will be discussed 
since they represent extensions of some of the previously reported analytic 
solutions. 

4.2 Analytic Solutions 

4.2.1 Shallow-angle Hypersonic TAV Climb 

Tauber and Adelmanl3 studied the flight mechanics of a shallow-angle 
hypersonic climb of a TAV on the way to orbit. The assumptions of the analysis 
were l)Y«lso that si 1 ; 2) constant acceleration ^ and dynamic 

pressure and 3) planar flight (V= constant). Under these conditions. 

Equations 1,2 reduce to o) 

® V “ *> (a- y) + »«*c«*eoj *f> * 
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The Tauber and Adelman analysis neglected Coriolis tenns so that Equations 15, 
16 become dlV _ T _ ^ 
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The study Went on to determine several analytic results including the following 
expression for fuel fraction as a function of propulsion average specific 
impulse and vehicle average lift/drag ratio: 

In the present study, a corresponding approximate result has been developed 
where the Coriolis effect is included. Consider the following equivalent 
assumptions to those made in the Tauber and Adelman studyrl) V =constant<*l and 
cos'® *1 ; 2) planar flight c = 0 and cojftant heading Y . In this case, Equations 
15, 16 including the Coriolis effect become 
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As in the previous study, the resonable assumption is made of average values of 
L/D and specific inpulse during hypersonic ascent. Further, it is here assumed 
that cos <P& constant, which is a restrictive assumption implying either small 
latitude changes or flight near the equator. With these assumptions (all but 
the assumption on cos^ are also assumptions in the previous study), Equation 21 
can be integrated with the result 
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This result is the extension of the earlier result (Equation 17) where an 
approximation to the Coriolis effect is now included. It should be noted that 
the approximation in Equation 21 represents a maximum Coriolis effect^ Figure 1 
shows plots for Equations 17 and 2A, and it indicates that the Coriolis effect 
reduces the average specific impulse requirement for fixed mf/m^ b y 4 to 

7%. For flight at higher latitudes (or greater latitude changes), the reduction 
will be less pronounced, i.e., the corresponding plot v/ill lie between the two 
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4.2.2 Hypersonic Gliding Turn 


A second extension to previous results was achieved in terns of including 
Coriolis effects (in an approximate way) in a basic relation between heading 
change and velocity change during a hypersonic gliding turn. Slyel4 developed 


the relation 


<S> 


where heading 0 and velocity V = Vi at the beginning of the turn, D = drag 
force, y = side (lateral force) = (L/D)sin«* if produced by vehicle bank angle «: 
Slye derived this result under the conditions of no thrust, shallow flight path 
angle (K«l), mgsinTr«D, constant L/D, and constant Y/D. This relation was 
used by Tauber and Yang25 to derive an equilibrium glide relation between 
altitude and heading, and then to produce a solution for constant ban); angle 
which minimizes velocity loss in the turn. 

The result given by Equation 23 is obtained by neglecting two terms in the 
complete lateral-force differential equation of motion (Equation 3). The first 
neglected term is the curved flight path centrifugal term ( )!2 coXcns V* ) > 

and neglecting it amounts to assuming negligible magnitude, perhaps due to 
restriction to near equatorial locations (small CP ). The second neglected term 
is the Coriolis term [ ao/V Ck>«Ycos<tsi»V — <t> ]. The present study obtains an 
approximate result containing the Coriolis term. Under the conditions assumed 
by Slye, Equations 1-3 became 
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which, upon substitution into Equations 3^- A, gives 
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Dividing Equation 29 by Equation 27 gives 

~ 'y 2i£c«>s0»n$ •*• aJ=^*L v (KcoiqiiA v-s>«c& 3 

This equation may be integrated to give very approximate results if average 
values (over the integration variable velocity V) are assumed for the 
products/differences of trigonometric functions of heading if> and latitude^. 

The results ^ ^ u (Vjf (jrcrts* V-s^ct] 

Further, an order-of-magnitude analysis shows that in 

equation 30, the Coriolis term is second order, while the (v2 - v? ) term is 

T 

first-order, comparable to the ln(V$/\{-) term. Thus, Equation 30 simplifies to 
(jj£> <f£. <f c — 

Figure 2 shows comparable plots of Equations 23 and 31 where average latitude, 
heading values of 30° and 20°, respectively, have been assumed in Equation 
31. Unfortunately, the subsequent relationships derived from Equation 23 by 
Tauber and Yang25 are difficult to parallel when beginning with Equation 31. 

Ikawa36 has examined rotating-Earth effects on AOTV trajectory 
simulations, and has concluded that they are signiificant (see conclusions of 
the present study below) . 


5. Conclusions and Recommendations 

The research performed during this effort involved both an evaluation and 
organizing of a rather wide literature and the techniques and solutions 
presented therein, and attempts to develop new analytical solution methods for 
parametric analysis of vehicle performance and design. The literature 
evaluation was very successful, but only limited success was achieved with 
extending previous analytical solution methods, and no success with developing 
completely new methods. 

It is believed the most important contribution of this research is the 
following set of conclusion: 

1) The best formulation of the governing equations for current missions of 
interest appears to be that of Equations 12 - 14: the plane of 
reference is the initial heading plane, the motion is Earth-relative, 



and the motion of the maneuvering hypersonic vehicle is easily 
expressed and visualized in crossrange/downrange coordinates. 

2) These equations can and should be implemented (and optimized for solution 

speed) on modern, high-speed, high-resolution graphics engineering 
workstations. 

3) Parametric design requires high-speed, fast-turnaround, user-interactive 

workstation-hosted implementation of these equations for trajectory 
simulations. This appears to be the only way to conveniently proceed. 
— £4) inclusion of Coriolis effects by the present study along with the 
study of Ikawa36 indicates previous analytic results could have 
significant errors in design and performance concllusions, due to 
neglect of these effects. 
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